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Abstract 

Using the gradient expansion method of Rigopoulos, Shellard and van 
Tent which treats cosmological perturbations as gradients on top of a 
homogeneous and isotropic FRW background, we study the production 
of nongaussianities in the roulette model of inflation. Investigating a 
number of trajectories within this two-field model of inflation, we find 
that while the superhorizon infiuence of the isocurvature modes on the 
curvature bispectrum produces nonzero contribution to /atl, the effect is 
negligible next to the standard infiationary prediction |/]vl| ~ Us — 1. 
This is the case in both the squeezed and equilateral configurations of the 
bispectrum, although the former is slightly larger in the trajectories under 
consideration. 



1 Introduction 

Recent advances in vacuum stabilization in string theory, most notably the 
KKLT [T] and the large volume [5] varieties, have put stringy realizations of 
cosmic inflation on a much firmer footing. Only when all moduli have been 
stabilized does it make sense to study which of them might be candidates for 
inflation. String theory provides an abundance of moduli fields, for example the 
size and shape of the extra dimensions, and thus many potential opportunities 
for inflation. These include single-field models, such as ref. [3j, and multiple- 
field models such as the racetrack model [3], the N-flation scenario [5J, or the 
swiss-cheese scenario [6]. Multiple field inflation has the advantage of more 
easily providing some of the nonstandard observational features that one would 
like for helping to discriminate between theories. The extra degrees of freedom 
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allow the production of isocurvature modes (perturbations transverse to the 
classical trajectory) and it has become apparent [7] that these may give rise to 
large nongaussianities in the cosmic microwave background (CMB) temperature 
fluctuations. 

One such model is the roulette scenario 0, in which the Calabi-Yau mani- 
fold (the compactified extra dimensions of string theory) relaxes from an initial 
excited state towards a minimum of its potential. The large volume compact- 
ification that is used ensures that this minimum exists for large ranges of the 
microscopic parameters. In addition, unlike KKLT, it does not require tuning 
the constant term in the superpotential to very small values, and it gives a 
natural expansion parameter, the inverse volume 1/V, providing a controlled 
a' expansion. In this particular model, the last four-cycle and corresponding 
axionic partner to relax act as slow rolling scalar fields which drive the final 
stage of inflation. 

Specializing to a specific set of microscopic parameters, we explored the 
various trajectories of inflation available within the context of this two-field 
model. The prediction of observable parameter values were in accordance with 
known results from CMB and large scale structure survey data. The influence 
of isocurvature perturbations, which seed inhomogeneieties between the species 
of fields driving inflation, was furthermore found to be quite important. 

In addition, recent claims of detection of nongaussianity in the WMAP CMB 
data, specifically a nonvanishing nonlinearity parameter fiqL [H] have sparked 
a large interest in deviations from gaussianity in the spectrum of primordial 
fluctuations as an additional observable that must be predicted by a successful 
theory of the early universe. Given that precision measurements of Jnl from 
experiments such as Planck will soon be available, it is all the more impor- 
tant that the mechanisms governing the production and evolution of primordial 
nongaussianities be well understood. 

We will first provide an overview of the roulette model, followed by a discus- 
sion of primordial nongaussianities from inflation. In Section [5] we will give the 
main results of our paper, which are the numerical caluclations of nongaussian- 
ities in the curvature perturbation. The important aspects of our calculations, 
which follow [TU], are presented in Appendix A. 

2 Roulette (Kahler Moduli) inflation 

The roulette model is a string theoretic inflationary scenario set in the context 
of a Type IIB large volume compactification. Although there may be evolution 
of several Kahler moduli, the observable part of inflation is governed by last 
(and lightest) one to relax. Since the earlier-evolving moduli stabilize to deep 
minuma, they rapidly decouple from the dynamics [8], before the final 60 e- 
foldings. The name "roulette" comes from the cyclic shape of the potential, 
resembling a roulette table whose grooves are the minima toward which the 
inflaton eventually relaxes (Figure [l}. During inflation, the F-term potential 
of the large volume compactification is flat enough to allow slow-rolling over 
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sizeable patches of field space. Reheating, which we will not address here, occurs 
when the inflaton fields oscillate at the bottom of the potential. This model was 
first proposed as a single- field inflation model by Conlon and Quevedo in [3], 
and subsequently generalized to include the axion as a second inflaton field by 
Bond, Kofman, Prokushkin and Vaudrevange in ref. [Sj. As in ref. [5], we use 
the large- volume compactification in which the 10 spacetime dimensions of 
type IIB string theory are separated into a 4-dimensional noncompact spacetime 
and a conformally Calabi-Yau 3-fold. 




Figure 1: The potential ([T]) as of function of the volume modulus r and its 
axionic partner 9 for parameters in Table 2. The potential is periodic in the 6 
direction. 

After minimizing the F-term potential in the large volume compactification 
with respect to the axio-dilaton and all Kahler moduli but T2 = T2 + i92, the 
potential can be expressed as a function of r = r2 and 9 = 62- After expanding 
in powers of 1/V this potential reduces to jHj: 

V ^ ^ m,a2A2re-^^^ cos(a,g) ^ ^ 

3aA2V 

where ISV is the uplifting contribution, of order \/V^ ^3], adjusted so that 
that y = at its minimum. It is crucial for the naturalness of Kahler moduli 
inflation that Ay depends only very weakly upon T2 in the case of interest, 
where ri ^ T2. Because of this, AV^ is nearly constant during inflation, and 
the slow-roll condition on T2 is easily satisfied. Fig. [T] shows the form of the 
potential. 
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For notational purposes we define cf)^ = t = T2 and (j)"^ = 6 = 62- The kinetic 
term takes the diagonal form: 

Ckin = \K22SABd^cl>^d^cl)'', (2) 

with the (2, 2) component of the Kahler metric given by [5] 

_ 3aA2[2V + e + 6aA2Tf'] 



In the original Kahler moduli inflation model 3 the Standard Model was 
confined to a D7-brane that wraps the inflaton cycle. It has recently been 
pointed out, however, 11] that string loop corrections contribute to the Kahler 
potential in such a way that the exponential flatness of V is destroyed. This 
would give large contributions to the rj parameter, thereby preventing slow roll 
from occuring. Ref. [11 does point out that this can be remedied by removing 
the D7 brane from the inflating cycle, which has the disadvantage of complicat- 
ing the reheating process. This could proceed via a mechanism akin to the ones 
studied in multiple-throat inflationary models in warped compactifications. See 
for instance ref. |12j . Given that we are only concerned with the inflationary 
phase in the present Paper, we will take the potential ^ as is. 



3 Inflation and nongaussianities 

Slow-roll inflationary scenarios generically predict a near-scale invariant spec- 
trum of primordial perturbations, with nearly Gaussian statistics. This is ex- 
pected, since individual quantum fluctuations can be treated as independent 
results and should therefore be Gaussian by virtue of the central limit theorem. 
However, more complex field interactions, along with interactions with gravity, 
which is inherently non-linear, are sure to produce at least small deviations from 
Gaussianity. Due to the stochastic nature of these perturbations, it is therefore 
natural to develop statistical tools to compare predictions of the theory with 
CMB observations. Wick's theorem tells us that the even moments (2n-point 
correlators) of a linear field or distribution 0^ decompose into a sum over the 
permutations of two-point correlators, whereas the odd moments vanish. The 
measurement of deviations from gaussianity in a field (/)(x), can therefore be 
made through the bispectrum (the Fourier transform of the three-point corre- 
lator), and through the connected part of the trispectrum, that is, the part of 
the Fourier transform of the four-point correlator that cannot be decomposed 
into products of the power spectrum. 

It is common to parameterize the small deviations from gaussianity in terms 
of their effect on the Bardeen potential <f> through /jv l [H] 

<i>(x) $^(x) + f^^ ($i(x) - ($i(x))) , (4) 

where $i is a purely gaussian random field with (^l) = 0. Even for large /jvl 
this parameterisation is sufficient, given that the fiuctuations $ are order ~ 
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Table 1: Some recent 95% CL estimates of Jnl using WMAP data. 



Komatsu et al. (WMAP 1-year) 


-58 < Jml < 134 


Creminelli et al. (WMAP 3-year) ^ 


-36 < fNL < 100 


Yadav and Wandelt (WMAP 3-year) ^ 


27 < fNL < 147 


Komatsu et al. (WMAP 5-year) [15] 


-9 < fNL < 111 



10^^. Already with the COBE observations it was shown that the nongaussian 
fraction of $ must be less than a few percent, fNLi^t)^^^ < 0.04 (see for 
example [13'). Subsequent measurements have tightened this limit to the level 
of fNLi^'i}^^^ < 0.003. Some recent published limits on fNL are shown in 
Table □ 

Some of these results [9] suggest that the CMB anisotropics exhibit measur- 
able deviations from Gaussian statistics. Whether or not these detections are 
confirmed, future observations such as the 9-year WMAP data and the Planck 
satellite data will provide stringent bounds on the primordial bispectrum, yield- 
ing additional parameters that any successful model of the early universe will 
have to match. Although nongaussianities may be measured from the contri- 
bution of any n > 2 connected n-point function, the 3-point correlator is the 
easiest to detect due to the smallness of the anisotropics. 

In the original ("local") ansatz fNL was taken to be a number, but 
by relating it to the bispectrum one sees that more generally it could be a 
function of the momenta fc^. Taking the Fourier transform and writing $(fc) = 
^L{k) + ^NL{k), is easy to show (noting that ($|,) vanishes identically) that 
the lowest order nonvanishing component of fNL may be written in terms of 
the bispectrum and power spectrum^ 

f I V , 1, , (^L(fcl)$L(fc2)^jVL(fc3)) 

Due to the delta function, the wave vectors form a triangle, and the k- 
dependence of eq. ^ can be expressed in terms of two ratios of momenta, 
for example fc2/fci and k^/ki, and an overall scale. Different mechanisms pro- 
duce bispectra that peak for differently shaped triangles; for example equilateral 
(fci = fc2 = fcs), or squeezed (fci <C ^2 = fca). The latter corresponds to the pre- 
diction of local ansatz. Single- field slow-roll inflation predicts nongaussianity 
of the local type [18], given that the dominant contribution to the bispectrum 
should come from the superhorizon influence of small k modes which act to 
"rescale" modes as they evolve toward the end of inflation. A rigorous expan- 
sion of the action to third order in perturbation theory is given in ref. jl8| . 



^ One must be careful with the sign of /jvl , which has been a source of some confusion in 
the hterature, due to the sign difference between the Bardeen potential and the gravitational 
potential (see appendix A2 of ref. [25].) We use the WMAP convention that positive Jml 
corresponds to positive bispectrum of 4?. 
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Other models such as ghost inflation and DBI inflation predict large Jnl for 
the equilateral configuration [19j . in which non-gaussianities are created before 
horizon-crossing. 

Our focus will be on tracking the perturbative curvature modes from the 
time they expand beyond the Hubble radius H~^, until the end of inflation. 
We will use the gradient expansion approach developed by Rigopoulos, Shellard 
and van Tent [20l [71 [10] . The assumptions of homogeneity and isotropy of the 
background FRW inflationary universe allow the use of the "long-wavelength" 
approximation, in which the gradient terms of the equations of motion may be 
dropped in the classical (unperturbed) equations of motion [TB] : 

An^ + 3A^ijn^ -NG^^Vb, (6) 

dtH = -^NHaU^, (7) 

where the Hubble rate is 

H'' = ^ (/Ikin + V) , (8) 

on scales larger than the Hubble length H~^. As a consequence, the inclusion 
of metric and field perturbations simply amounts to the inclusion of gradient 
terms on top of the background, whose equations of motion are computed from 
the full field equations of motion. The advantage of this approach is twofold: 
no slow-roll approximation is needed to find and solve the equations of motion, 
and the resulting equations are exact (non-perturbative) results. Quantitative 
computation of power spectra and bispectra was done via a perturbative ex- 
pansion of these gradient equations of motion. The important results of [TOj 
that we made use of are presented in Appendix A. Here we will only sketch the 
method and refer the reader interested in the details to the Appendix. 

For the two-field case such as we consider here, the dynamical degrees of 
freedom for the metric fiuctuations can be reduced to (C^, C^, 0^) where is the 
adiabatic curvature fiuctuation, is the isocurvature fluctuation, and 9^ = C2- 
{9^ is like the canonical momentum conjugate to (^^, and the evolution equations 
are first order in derivatives, similar to Hamilton's equations of motion. 9^ is not 
an independent degree of freedom, but is constrained in terms of the others.) 
The formalism solves for the gradient of the three fields, Via = di{C} ^ ^ 9"^) , 
where a = 1,2,3. These quantities can furthermore be expanded order- by-order 
in the cosmological perturbation, Via = + v^-^ -I- . . . . Higher order terms are 
sourced by the next lowest order ones, through a master equation of the form 

Via{t, X) + Aabit., X)Vib{t, x) = 0. (9) 

where the matrix Aab (see (IA.18|) is determined by various slow roll parameters. 
The lowest order source term from which higher order fiuctuations follow is 
deduced by the method of stochastic quantization. This information is encoded 
in a matrix xj^^^, eq. (|A.20|) . The solutions of eq. Q for vH'' and v^"^^ which 
we will need for computing the spectrum and bispectrum can be expressed in 
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terms of a Green's function Gab(t — t') which is the sohition to 

-^Gabit, t') + 4°) {t)G,,{t, t') = S{t ~ t'). (10) 

where is the matrix Aac evaluated using just the homogeneous background 
solution. The main technical difficulty then is in computing the Green's func- 
tion. The fluctuation can then be computed through wim(fc, t)Gofc(t, t* + 
lnc)X^^(fc,i» + Inc) (see eq. (jA.31|) ) and similarly vlf* is given by eq. (|A.37|) . 

The nongaussianity parameter /jvl is determined by and v^'^'^ through eqs. 
(|A.35|) and (|A.36|) of the appendix. 

For the roulette model, we must solve these equations numerically. Analytic 
results have been developed {e.g., in ref. [7]) in the context of this formalism, 
but only within certain constrained limits. 

4 Numerical method and results 

In this section we will consider a number of different parameter sets for the 
model, indicated in Table 2. Special attention will be given to the first of these, 
for which we illustrate the different possibilities depending on which inflationary 
trajectory is followed {i.e., the dependence on the initial conditions). In each 
case, we evolved the light fields (r, 9) starting from rest, until the end of infiation, 
which we took to be the point at which the slow roll parameter e grew to e = 1 . 
For each trajectory, we solved for the perturbation amplitudes v^^^^ to first and 
second order using the gradient method described above. The power spectrum, 
scalar spectral index and nonlinearity parameter /jvl thus calculated are pre- 
sented below for a variety of inflationary trajectories and parameter choices. 
We also discuss the superhorizon influence of the isocurvature modes on these 
quantities in what follows. 

In contrast with the KKLT compactifiation, the large volume scenario places 
no strong restrictions on the value of Wq in the effective field theory [5]. One 
therefore has considerable freedom to vary the parameters of the potential ([!]). 
We first verified the results of ref. [8] as a check on methodology. We solved 



Table 2: Parameter sets used for numerical simulations of the roulette model. 
The first corresponds to parameter set 1 of [8], but with V adjusted to meet 
COBE normalization and was used in all simulations unless otherwise indicated. 



Set 


Wo 


02 


A2 


A2 


a 


e 


V 


1 


300 


27r/3 


0.1 


1 


1/9 V2 


0.5 


8 X 10* 


2 


6 X 10"* 


27r/30 


0.1 


1 


1/9 V2 


0.5 


10* 


3 


4 X 10^ 


tt/IOO 


0.1 


1 


1/9 V2 


0.5 


109 


4 


200 




0.1 


1 


1/9 V2 


0.5 


10^ 



7 



for the inflation trajectories starting from a variety of initial conditions and 
verified that slow roll was generically obtained, as found in section 6 of [8]. 
Some examples are illustrated in Figure [51 Solution A, the "r- valley" trajectory 
of Conlon and Quevedo [3], effectively corresponds to single-field inflation, as 
the fields start with already minimized. Figure [3] shows a more detailed plot 
of one of the trajectories, superimposed on a contour plot of the potential. 
Values of e were consistently very small, with loge ^ —13 at the COBE scale. 
Typical values of the tensor-to-scalar ratio produced by the fields were therefore 
r ~ 3.5 X 10^1^ 

We focused attention on parameter set 1 of [5], with the modification that 
the volume V was tuned to achieve COBE normalisation (V = 8 x lO^^f; see 
Table[2]), i-e., the power spectrum is Ps ~ 4 x lO^^*' on COBE scales, in order 
to have a realistic example. This corresponds to an inflationary energy scale 
yi/4 ^ iQi3 QeV, giving a duration of 52-55 observable e-foldings of inflation, 
assuming reheating temperatures of 10"'^°-10^^ GeV. Although a more generic 
method of normalization is to rescale the potential by an overall factor, the 
dependence on 1/V of both terms in the potential gives a way to adjust its 
magnitude without introducing additional parameters. 



50 r 




Figure 2: Various field trajectories for different initial field configurations (red 
dots). The numbers beside each curve are the number of e-foldings before slow 
roll breaks down and e exceedes 1. The potential here used the parameters from 
Table [5] The labeled trajectories A through G correspond to those listed in 
Tabled 
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Figure 3: Detailed plot of trajectory C from Figure[2l superimposed on a contour 
plot of the inflaton potential. Red dots represent time increments in e-foldings. 



4.1 Isocurvature Perturbations 

In the case of more highly curved trajectories, the isocurvature modes were 
found to have a large, positive effect on the power spectrum of adiabatic per- 
turbations, and consequently on the scalar spectral index n^, as recently shown 
in [21] . Large curvature in field space during the course of inflation resulted 
in a "projection" of the isocurvature modes onto the adiabatic direction. For 
some trajectories with long periods of curving, over 90% of the power spectrum 
originated from the isocurvature modes. Table [3] gives the proportion of the 
observable curvature power spectrum at the end of inflation that results from 
the influence of the isocurvature mode: 



(11) 

t,=55 



where the subscript = 55 indicates that these quantities were evaluated at 
the COBE scale, which wc take to be the modes that crossed the horizon 55 
e-foldings before the end of inflation. Ps.f . is the power spectrum computed with 
the effective single-field result 

whereas Pcxact was computed using (jA.34p . A second result of the influence of 
isocurvature modes (also discussed in ref. |2j) is a lower scalar spectral index 
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Table 3: Cosmological observables computed from chosen trajectories (see fig- 
ure [2]). The top part uses parameter set 1, whereas the bottom part displays 
results from the other sets (see Table [2]). /^^''' and f^^'^' are the nonlinearity 
parameter in the equilateral and squeezed configurations, respectively; N is the 
total number of e-foldings of inflation; n|-^- = 27y — 6e is the standard single-field 
scalar spectral index, whereas n™'*'^' is the exact index. p'^° is the proportion of 
the curvature power spectrum originating from the isocurvature modes during 
superhorizon evolution, calculated using pT|) . All observables are computed for 
the mode k that left the horizon 55 e-foldings before the end of inflation, and 
are evaluated at the end of inflation. 



Trajectory 


a2T 


Si20 


oequil. 


f-sqzd. 


N 




j^exact 


piso 


A 


48.4 


TT 








86 


0.965 


0.965 





B 


47.6 


7r/2 


0.0052 


0.0069 


62 


0.992 


0.976 


0.70 


C 


47.0 


7r/4 


0.0105 


0.0132 


72 


1.02 


0.940 


0.81 


D 


46.3 


7r/8 


0.0242 


0.0272 


77 


1.02 


0.920 


0.91 


E 


38.95 


0.36 


0.0688 


0.0698 


77 


1.04 


0.930 


0.90 


F 


38.93 


0.40 


0.0068 


0.0089 


103 


1.002 


0.948 


0.72 


G 


38 


0.38 


-0.0060 


-0.0060 


90 


0.964 


0.965 


0.02 


Param. set 


a2T 




oequil. 


f-sqzd. 


N 


s.f. 


j^exact 


piso 


2 


29.5 


7r/32 


0.0301 


0.0338 


77 


1.055 


0.909 


0.93 


3 


28.7 


7r/32 


0.0553 


0.0590 


69 


1.1045 


0.849 


0.97 


4 


33.2 


7r/32 


0.0404 


0.0442 


73 


1.074 


0.891 


0.93 


4 


35.0 


7r/8 


0.0060 


0.0082 


94 


1.002 


0.946 


0.73 



than would be naively expected from the single-field result = 2ry — 6e 

(= — 4e — 677" in the notation of Appendix A). The full scalar spectral index Ug 
at COBE scales was computed by taking the derivative of a cubic fit of In Pg , 
with the power spectrum evaluated from (|A.34p : 

(exact) _ C^lH-Pcxact _ In Poxact 

~ d\nk " di 

Here we have used the fact that in the gauge NH = 1 with H approximately 
constant, d/dhik = d/d\n{aH)^ ~ H~^d/dt. Table [3] gives a comparison of 
both methods of computing Hs- Our results indicate that a significantly larger 
power spectrum, along with a generically red-tilted spectrum is an expected 
result of curved trajectories in roulette inflation. This is of particular interest, 
given that the most recent cosmological data favor a scalar spectral index of 
rig = 0.96 [Tg. 



(13) 
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10 20 30 40 50 60 70 80 

t 

Figure 4: Green's functions Gab{t,t') for t' = 60 e-folds before the end of infla- 
tion, for roulette inflation with the parameter set 1 and trajectory C. Only the 
nonzero components are labeled. 

4.2 Nongaussianities in the roulette model 

We now turn to the subject of nongaussianity in the roulette model. Because it 
gives examples of highly curved field trajectories, one might have hoped to find 
observably large levels of nongaussianity coming from the isocurvature modes. 
However, the numerics do not bear out this expectation, as we now describe. 

We tested the algorithm for computing f^^ from superhorizon evolution of 
perturbation modes on the two-field quadratic inflation model considered in ref. 
|10| . verified its results. We then analyzed the nonlinear mode evolution for a 
variety of roulette infiation trajectories, for modes k corresponding to a range 
of horizon exit times = Ink/H^ before the end of inflation. 

The Green's function in eq. (|10p was found by solving the ODE numerically 
in matrix form, as a function of t. This was done once per time step t', giving 
a 3 X 3 X Af X Af dimensional array, where M corresponds to the number of 
discrete time steps sampled in the evolution (typically around 1000). Figure |4] 
shows this behavior for trajectory C of flgure 2. 

Inflation in the r valley, corresponding to the effective one-field scenario of 
Conlon and Quevedo [3 produced no nongaussianities originating from super- 
horizon interaction between scalar modes {/nl ^ 10^^°, where the deviation 
from zero can be attributed to numerical uncertainty). This is not surprising, 
since it is the coupling between curvature and isocurvature modes that is ex- 
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pected to generate large bispectraH In more complex inflationary trajectories 
with sufficient curving in field space, however, we found that values of /^vl be- 
tween —0.01 and 0.02 were quite generically produced. But we did not find any 
examples which produced /nl outside of the small range 

l/iVLl < 0.1 (14) 

This is far below the level of sensitivity foreseen by the PLANCK experiment, 
for example. 

Figure [5] illustrates the time-dependence of /nl at a series of wave- numbers 
k from horizon exit to the end of inflation, for the representative trajectory C 
in Figure [H), while figure [6] shows the /c-dependence of Jnl after it has stopped 
evolving. The behaviour of /nl shown is typical for trajectories that continued 
to curve during most of the period of observable inflation: they produce slightly 
more pronounced values of /nl during the curved part of the motion, but 
these values quickly descend to ~ 10~^ by the end of inflation. The two-field 
model studied in jlO] also exhibits this behaviour: /jvl descends to zero as 
the trajectory straightens out at the end of infiation, resulting in a bispectrum 
below the level of measurable sensitivity. Table [3] gives some computed values 
of /^/^''' (where the ki form an equilateral triangle) at scales k that crossed the 
horizon at t^{k) ~ 55 e- foldings, as well as for squeezed triangles, with fci, 

k2 and k^ corresponding to i»(fci) = 60 and (^2) = t^lk^) = 55 respectively. 
Trajectories labeled A through G are shown in figure [21 The bottom part of 
table [3] shows results for different parameter sets (given in Table [5]) . Values of 
/jVL are taken at the end of inflation. 

Features in fNhik) can be understood as being due to curvature of the tra- 
jectory at the time of horizon exit of particular modes. Modes that experienced 
more curving after horizon exit (i.e., those which exited the horizon earlier) pro- 
duced larger magnitude /jvl than those which experienced no curving. These 
features in the bispectrum are completely absent in single-field inflation. One 
feature that was common to curved roulette trajectories was a slightly larger 
bispectrum on large scales, due to the modes which left the horizon before turn- 
ing in field space occurred. This occurs because only superhorizon-scale modes 
can experience growth due to the coupling between isocurvature and adiabatic 
modes (this effect can also be seen in Table [31 which shows the correlation be- 
tween large isocurvature contributions to the power and larger values of /nl)- 
Smaller scales which exit the horizon later undergo less such growth. Thus a 
simultaneous detection of larger /jv l in the CMB and smaller primordial non- 
linearity in large-scale structure may be a way to detect this type of result, if 
its magnitude can be enhanced to an observable level. 

One somewhat nonstandard feature of the roulette model, relative to simpler 
two-field models, is the nontrivial field metric K22 (eq. ([3]) which multiplies the 
kinetic term of the axionic direction 9. We found that it had no substantial 
effect on the shape, size, or magnitude of the bispectrum. To test this, we tried 

■^Recall that contributions to /jvl from the adiabatic fluctuations should give values of 
order (n^ - 1) [18]. 



12 
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Figure 5: Time evolution of the non-linearity parameter /tvl for various modes, 
for trajectory C in the equilateral configuration. The beginning of each curve 
corresponds to t ^ t^,, the time of horizon crossing of the mode. Here, the x-axis 
is time until the end of inflation, in e-folds. 

replacing K22 by a constant during inflation, and found that the spectrum of 
observables remained roughly unchanged. 

As illustrated in Table [3l the shape-dependence of the observed nongaus- 
sianities is as expected from this type of model [19j , since nongaussianities from 
superhorizon evolution are larger in the squeezed configuration, whereas those 
generated on subhorizon scales (like in DBI inflation) are larger in the equi- 
lateral configuration. In all cases, we find that /tvl is larger for the squeezed 
configuration. 

5 Discussion 

We have studied a model of Kahler moduli inflation built from a realistic con- 
struction of Type IIB string theory, using the formalism developed in [20l [71 [10] . 
We confirm previous results for the power spectrum and the superhorizon influ- 
ence of isocurvature modes for typical inflationary trajectories. The main new 
undertaking is a search for nongaussian perturbations from highly curved tra- 
jectories. Such deviations from gaussianity can originate from the superhorizon 
evolution of the second-order curvature perturbation and its interaction with 
isocurvature perturbations. 

When the full spectrum is considered, roulette inflation predicts a smooth 
power spectrum with a slight red-tilt, in excellent agreement with estimates 
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Figure 6: Wave- number dependence of Jnl for trajectory C. The x-axis is the 
time of horizon crossing in e-folds of the corresponding mode, where i» = is 
the end of inflation. In the case of f^^*^' we took ii* = 70 e-foldings. This is 
the reason that = 70) = fj^L^'it* = 70). 
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based on the latest WMAP data. In addition, as shown by [21], the superhori- 
zon influence of isocurvature modes can come to dominate the scalar curvature 
power spectrum via the relation (jA.15|) . This is the case for inflationary tra- 
jectories with large curvature, which can be pictured as a "projection" of the 
isocurvature perturbation modes onto the adiabatic direction. This contribu- 
tion, which can account for over 90% of the curvature power spectrum, must 
therefore be considered in the context of these multiple-field inflationary models. 

Concerning the issue of nongaussianity, we did not find examples in which 
fNL^ produced by the mechanism considered in section 21 could be large enough 
to be observable by future missions; we obtained results which did not exceed 
the level of f^L expected in conventional single-field inflation models. It was 
hoped that trajectories with very sharp curving in field space would have yielded 
larger /jvl values, but two considerations could make this difficult. First, if the 
trajectory relaxes to a straight, effectively single-field form before the end of 
inflation, Jml damps to zero, as illustrated in figure 1 of ref. [101 ■ We found 
trajectories with moderate curving until the end of inflation, which prevented 
the total erasure of Jnl- Multiple-fleld ( "multi-brid" ) hybrid inflation models 
such as discussed in refs. |22l [23] use this effect to produce sizable nongaussian- 
ities. Second, large deviation from Gaussianity from this method is correlated 
with large isocurvature modes being projected onto the curvature direction. If 
one is able to find models with large /nl, it must be verified that the power- 
spectrum is not overly amplified and distorted by this effect. This phenomenon 
is also seen to play an important role in nongaussianity generated by the wa- 
terfall fields in hybrid infiation models [21] . 

However, it is possible that a more complete search of the model's parameter 
space would reveal examples with larger values of /nl- In that case, running 
bispectra {i.e., which depend on the magnitude of the average scale of fci, as 
opposed to the shape-dependence), such as we find in the roulette model, could 
provide an interesting discriminator between models. Recent developments in 
the detection of nongaussianity via large-scale structure, which would probe 
Inl at smaller scales, promise to give additional observational handles on such 
a dependence [251 [53] . 

It may be interesting to examine other such moduli inflation scenarios that 
arise once the assumptions of a strict hierarchy of scale are relaxed (for instance, 
the model of ref. in which the second dynamical field is the inverse overall 
volume, rather than the axionic partner). A more general study of two-field 
inflationary models with such an exponentially flat potential would furthermore 
reveal how generic the above-mentioned behavior is. 
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A Second order quantities in the gradient ex- 
pansion formalism 

We present here the important resuhs of [10] , which were used in the computa- 
tion of the first- and second-order quantities that formed our main results. 

Regardless of the specific inflationary model, we may write the action for 
a general multifield potential V{(j)A) {A = 1,2,...) with minimal gravitational 
couplings as 

j (^^R - ^GABd^.c/y^d'^q^'' - 1/(0^)) . (A.l) 

Here R is the Ricci scalar, Gab is the metric in field space and Mp — 1/ VSnG is 
the reduced Planck mass. We work in units where Mp = 1. At the homogeneous 
level, and before specifying a choice of spacelike slicing, the FRW metric is 
ds^ = —N'^{t)dt^ -\- a^{t)dx^, and varying with respect to the fields (fiA, a{t) 
and N{t) gives the equations of motion for the scalar fields and the Friedmann 
equations. Here N{t) is the time-lapse function, and a(t) is the scale factor, from 
which we define the Hubble parameter H{t) = aj {No). To simplify calculations 
during inflation, we make the coordinate choice 

< = ln(a), (A.2) 

so N{t) = H^^ ~ constant during inflation. Dotted fields will from now on 
represent differentiation with respect to this time parameter. 

Working within the long- wavelength approximation '20j, we assume that 
the fields are homogeneous and isotropic within the horizon. The kinetic term 
becomes: 

CMn = ^GAsH-^n^, (A.3) 

having defined the velocities 

= ^ = H^^ (A.4) 

We will furthermore define covariant differentiation of a field that transforms 
as a vector within field space [20] : 

DbL^ = dBL^ + TicL'' (A.5) 
DbLa = dsLA-T^Lc (A.6) 
D,,L^ = d^L^ + ricd^.4>''L^ (A.7) 

where Tg^j is the connection defined through the metric Gab- Henceforth, 
uppercase latin indices A,B,C,... will represent the various fields, greek will 
represent spacetime indices, and i,j,k will be spatial indices. Since we are 
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interested in a two-field inflation model, A,B,C,... = 1,2. In the roulette 
model there is only one independent connection coefficient: 

^rr-Ke-^er- ^66 - V + 3aAr3/2 ) (A.8j 

^le^Kr = (A.9) 

For the analysis of the power spectrum and bispectrum, it will furthermore 
be useful to define the orthonormal basis e^, where m — 1,2: 

4 = e^ = eABef (A.IO) 

where 11 = ^Ilyill"^ and eAs is the antisymmetric tensor, is tangent to the 
classical field trajectory, whereas is orthogonal. Note that lower and raised 
indices m, n are equivalent. 

The scalar field equations of motion are given in eqs. ([8][7]). We numerically 
integrated them to determine the infiationary trajectories. The formalism of 
PUI [71 [TU] which we follow makes extensive use of the "slow roll" parameters. 



kin 



7y" 



2i?2 ij2 

3gn^ + G^^dBV 

e^VA 
HU 



X = ^ + ^ + 

S.m = + 3(e - 77ll)(5ml - 377^5m2 

= a, (A.ii) 

where V22 and Vmi are defined in the orthonormal basis (jA.lOp . such that Vmn = 
e^e^V-AB, with the covariant derivatives over the field metric defined above. 
These quantities are nonlinear, depend on both t and x, and are not assumed 
to be small, although they are small in the slow-roll regime. It should be noted 
that some of them are unintuitively named; for example 77^ is proportional to 
the slope of the potential, in the direction orthogonal to the trajectory, rather 
than a curvature, and the relation of 77" to the usual slow-roll parameters, with 
respect to the adiabatic direction, is Ty" = —77 -j- e. Nevertheless we will keep 
this notation for ease of comparison with ref. [7] . The e parameter does agree 
with the conventional e (defined with respect to the slope along the adiabatic 
direction). 
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The particular choice of variables we will use to describe the fluctuations in 
this formalism are 

C^(t,x) = efl(t,x)ajna(t,x) - J. - d,cl>^{t,K) (A.12) 

^2e(t,x) 

which have the property of being invariant under long-wavelength changes of 
time-slicing {t,x) {i,x) [2D].This quantity can be projected onto the field 
basis (|A.10|) : 

Cnt, x) = S„,i d,\na- ^e„,A 9,0^. (A.13) 

These simplify in the gauge t — In a, where 9^ Ina = 0. At first order, Q is the 
spatial gradient of the usual curvature perturbation, whereas (^f corresponds to 
the isocurvature perturbation. The gradients are combined, along with their 
respective velocities 6*™ = 9fC™, into a 3-component vector, 

V^a^{Cl,Cf,Off, (A.14) 

The would-be fourth component is not independent, but is determined to be 

9l = 2rj^Cf (A.15) 

by the constraint equations [1^ which may be derived from the Einstein equa- 
tions and the definition of CI") noting that Dt{di4>^) = Di^NU-^): 

d,\nH = eC^ (A. 16) 

The relationship between 6l and Cf is nothing more than the well-known "con- 
servation" law of the curvature perturbation. This is valid to all orders, as 
shown in ref . [25 . 

Combining these with the equations of motion ©-([7]), the full nonlinear 
evolution equations may be written in the compact form: 

Via{t,x) + Aabit,x.)v,b{t,:ic) = 0. (A.17) 

The matrix A is a function of the parameters defined in eq. (|A.lip [TU] : 

/ -2ri^ \ 

A = 

\0 3x + 2e2+4er?ll +4(77^)2+^11 -2ei?2ii2 3 + e + 27/11 / 

Its dominant components are A33 = 3 and A23 = —1. The only explicit de- 
pendence on the curvature of the field manifold in A is the term — 2ei?2ii2 = 
—2ee2 ef Ci 62 Rabcd, but we found that this is negligible (~ 10^^) in roulette 
inflation. 
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The next step is to solve this system of equations perturbatively. Eq. (|A.17[) 

(n) 

can be expanded into a hierarchy of hnear perturbation equations for v^^ , each 
sourced by the previous order. Since we are interested in superhorizon evolu- 
tion, it is reasonable to take the first-order perturbations to be sourced by a 
linear perturbation , which encodes the effect of quantum fluctuations at 
short wavelengths providing the initial values for the long-wavelength modes of 
interest at horizon crossing. Refs. [3 [501 HH] show that the source term having 
the right properties is 

■^-^ W{k)Xi]lal{k)the^'^- + c.c, (A.19) 

where the creation operator has the standard commutator [am(k), a|j(k')] = 
Smn^^^^O^ — k'). Superscripts in parentheses indicate the expansion order in 
perturbation theory. 

The matrix of linear solutions around horizon crossing Xam is the slow-roll 
solution of ref . (lOj in which it is argued that deviations from linearity on sub- 
horizon scales should be slow-roll suppressed: 




^- = -7137^7= 1. (A.20) 



The factor l/v2e comes from the definition of C, and the amplitude H is the 
result we expect from perturbations at horizon crossing. 

The window function yV{t, k) is designed to source only the superhorizon 
modes, and the final results must be independent of its exact shape. It is 
convenient to use a Heaviside step function, >V(t, /c) = Q{kR — 1), that has 
support only on scales R = {c/aH) — {c/H)e~* (recall we are in the gauge 
t = lna) sufficiently larger than the Hubble radius, where c should be of order 
a few. Given that fluctuations that are generated on sub-horizon scales do 
not yet feel the effect of curvature, and therefore correspond to fluctuations 
in Minkowski space, it is reasonable to expect the spectrum of fluctuations on 
these scales to be Gaussian. 

Then 

m, k) = S{kR - 1) = (A.21) 

where is the time of horizon-crossing of mode k 

t^=\nk/H^. (A.22) 

Physical quantities are found to be independent of the exact value of c > 1. 
The first- and second-order equations can then be written: 

r;W(t,x)+4",)(t,x)z;(,^)(t,x) = hfj, (A.23) 
{t, x) + (^, {t, x) = - A« {t, x)«W {t, x). (A.24) 
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Here, A^^^^ — A^^i^^(t)d "^d^v^^ (note that 9 ^ is just multiplication by —k ^ in 
Fourier space.) A is given by [l0| : 

-6x - 2e?7ll - 2(7711 )^ - 2(r;^)2 

V -6 - 2,7" y 



/-2e2 _ 4ej7ll + 2(r;ll)2 - 2(77^)^ - 2,^" 
A32 -4677^ -2^^ 

V -2^^ JJ 

(A.25) 

where 

/ A321 \ 

A32 = -2ftei?2ii2 + A322 (A.26) 

V A323 / 

and 

A321 - -667711 -6(77^)2 -3ex (A.27) 
-4e3 - 10e277ll - 2e(77ll)2 

-6e(?7-L)2+8?7ll(?7-L)^-3e,e" 

A322 = -1267^-^-677117/-^ (A.28) 
+1277^X-6e2?7^ +4(?/-L)3 



^323 



-4ee^ - 277!!^^ + y 2("^2ii - 1^222) 

677^ - 267;^ + 477II77-L - 2^-^ (A.29) 



The first index in Aabc stands for the row, the second for the column in (|A.25|1 . 
and the third for the "depth" dimension of the array, represented here by a 
column vector for each Aat- Vimn is defined as Vimn = efe^e'^V-ABC- 

The above linear equations can then be solved with the aid of the Green's 
function which is the solution to the inhomogeneous equation 

^Ga6(t, t') + {t)Gcb{t, t') = S{t - t'). (A.30) 

with Gab{t,t) = Sab at equal times. This must be solved only once for each 
classical trajectory, which we do numerically on a grid in t, t' , a and b. Once 
Gabit,t') is known, the step-function form of W{t, k) simplifies the integration 
of the first order solution, 

^;«(fc,t) = /* dt'Gab{tX)MkX)xi2{k:t') 

J — 00 

= Gab{t,U+\Tic)X^ll{k,U+\nc), (A.31) 
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where we have defined the Fourier-space perturbation as: 

The second order sohition can be expressed, using the same method, as: 

via\^,t) ^-J dt'Gabit,t')Ai,cd{t')v'i^\^,t')d-'d^vfJ{^,t'). (A.33) 

To connect with observables one transforms the time coordinate t in the gauge 
of uniform expansion time sUces {NH = 1) to T{t,x) which describes uniform 
density sUces {dip = 0) so a{t) a{T,x). Then the curvature perturba- 
tion can be expressed as the total gradient of a scalar a, — i9ilna = dia, 
which allows observable scalar correlators to be expressed simply. Note that 
this result should be identical to results found using the SN formalism, given 
that the perturbation da in the uniform density gauge corresponds exactly to 
the perturbation in the number of e- folds 5N = 6\nd = ( [25] . 
The curvature power spectrum is: 

v{k,t) - ^{aa){k,t) = ^v['J^{k,t)v^^J^{k,t). (A.34) 

The scale-dependence of P comes as expected from the time-dependence of H in 
Xamik,t) (eq. (|A.20p ). which appears in v[]l{k,t) through eq. (|A.3ip . It should 
be stressed that the power spectrum here is complete and includes the effect of 
isocurvature perturbations. 

The leading contribution to the bispectrum comes from the expansion to 
second order in perturbation theory, 

(afe,afc,afc3)(2)(t) = {a^^^ a^^) ) (t) + {k^ ^ k,) + {k2 ^ k^) (A.35) 
= (2^)3,53(ki + k2 + kg) [/(fci, fca) + /(fc2, fcg) + /(fci, fcs)] 

where [10] 

f{k,k') ^ Qt;fl(fc,A:',i) +7y^««(fc,<)z;«(fc',i)) v^^2.{k,t)v[]^{k',t)+k ^ k' . 

(A.36) 

The second term in parentheses comes from the coordinate change t ^ T, and 
the first term is given by 

^Pl(fc,fc',i)^- /* At'G,a{t,t')Aa,S')v\^^{k,t')v'^^2M- (A.37) 



Numerically, we will find that this term dominates over the 77^W2m(^j ^)'^iri {k' , t) 
term in the roulette inflation model by five orders of magnitude. These are all 
the ingredients needed for evaluation of the nonlinearity parameter /atl [Z] > 

, (4?4>g) + (fcl-fc3) + (fc2^fc3) 

(a(i)a(i))fc, + {k, ^ k,) + (fe k^) ' ^ ' 
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